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ABSTRACT 

Bolometric corrections based on the optical-to-ultraviolet continuum spectrum of 
quasars are widely used to quantify their radiative output, although such estimates are 
affected by a myriad of uncertainties, such as the generally unknown line-of-sight angle 
to the central engine. In order to shed light on these issues, we investigate the state-of- 
the-art models of Hubeny et al. (2000) that describe the continuum spectrum of thin 
accretion discs and include relativistic effects. We explore the bolometric corrections as 
a function of mass accretion rates, black hole masses and viewing angles, restricted to 
the parameter space expected for type-1 quasars. We find that a nonlinear relationship 
log Lboi = A + B\og{XL\) with B < 0.9 is favoured by the models and becomes tighter 
as the wavelength decreases. We calculate from the model the bolometric corrections 
corresponding to the wavelengths A = 1450 A, 3000 A and 5100 A. In particular, for 
A = 3000 A we find A = 9.24 ± 0.77 and B = 0.81 ± 0.02. We demonstrate that the 
often-made assumption that quasars emit isotropically may lead to severe systematic 
errors in the determination of Lboi, when using the method of integrating the "big blue 
bump" spectrum. For a typical viewing angle of k, 30° to the quasar central engine, we 
obtain that the value of Lboi resulting from the isotropy assumption has a systematic 
error of « 30% high compared to the value of Lboi which incorporates the anisotropic 
emission of the accretion disc. These results are of direct relevance to observational 
determinations of the bolometric luminosities of quasars, and may be used to improve 
such estimates. 
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1 INTRODUCTION 

Quasars are important and interesting as the most lumi- 
nous non-exploding objects in the universe. They seem 
likely to exert a powerful influence on the evolution of 
the galaxies that host them, particularly through their ex- 
treme luminosities tha t can ionize and radiatively accelerate 
surrounding gas ( e.g., Crenshaw. Kraemer fc George _ 20031 : 
iBegelmanI |2004| : iDi Matteo. Springel fc HernquistI boOSi ). 



Accretion onto supermassive black holes powers quasars, 
and our ability to estimate the masses of these black holes 
has significantly advanced our understanding and revealed 
relationships among quasars and their h ost galaxies (e.g., 
iFerrarese fc Fordll2005l : IVestergaardll2009l ) . 

Other fundamental quasar properties include the bolo- 
metric luminosity and the accretion rate, often expressed as 
the Eddington fraction (the ratio of the bolometric luminos- 
ity to the Eddington luminosity). The ability to determine 
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the bolometric luminosity is often taken for granted, since 
the quasar continuum is observed directly in m ost cases and 
simp le scaling relationships are assumed (e.g., IVestergaardI 
l2009t ). How to make the bolometric corrections is not agreed 
upon, however, and the uncertainties are probably greater 
than those in determining black hole mass. Understanding 
how to properly make bolometric corrections and their un- 
certainties is important to future quasar investigations. 

The idea behind determining bolometric luminosity is 
simple: from one or more observations determine the total 
energy radiated by a quasar per second into its surround- 
ings. For objects like stars this is relatively easy, but that is 
not the case for quasars. First, quasars emit anisotropically, 
in a complicated manner, and the degree of this anisotropy 
depends on wavelength and likely varies significantly from 
object to object. Second, quasars emit at all wavebands, 
from radio through X-rays and beyond, many of which are 
difficult to observe. Third, portions of the spectrum repre- 
sent reprocessed photons; that is, some fraction of optical or 
UV photons can heat surrounding dust which in turn reradi- 
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ates at infrared wavelengths, so double counting can be an 
issue in theoretical calculations. Fourth, stars and stellar- 
heated dust can, in certain parts of the spectrum (e.g., far 
infrared, optical) contaminate quasar SEDs and be difficult 
to remove. 

These complications make it difficult to determine the 
bolometric luminosity of even extremely well-observed indi- 
vidual quasars, let alone that of a random quasar with only 
limited observations. Still, the bolometric luminosity is an 
important quantity and theoretical considerations can help 
us better understand the problem of bolometric corrections 
and may suggest improvements to the methods currently 
employed. 

One standard practice to determine bolometric correc- 
tions has been to use observations of quasar spectral energy 
distributions (SEDs) to determine the apparent bolometric 
luminosity, by integrating over all wavelengths, and then 
fitting relationships betwe en optical magnitudes and total 
bolometric luminosity (e.g.,|Elvis_et_al ] |l994l : lRichards et al.l 
I2OO6I ). It is generally assumed, using this approach, that the 
quasar is responsible for the light being emitted at all wave- 
lengths and that the emission is isotropic. The bolometric 
corrections derived generally show about an order of magni- 
tude of dispersion around the average value. This approach 
has generally only been used to fit simple linear relation- 
ships. 

Another standard practice (e.g., |Ne t zer et aL Il2003) has 
been to assume that the optical-ultraviolet portion of the 
quasar SED is energetically dominant and powers emission 
at other wavelengths. This approach avoids double counting, 
and is simpler in that only a portion of the SED, the so-called 
"big blue bump" thought to be emitted by a geometrically 
thin accretion disc, is important to know. Isotropic emission 
is generally assumed. This technique generates bolometric 
corrections not too dissimilar from those of the above tech- 
nique, and apparent uncertainties are also of about an order 
of magnitude. 

We can look toward accretion disc theory in an attempt 
to improve on these simple approaches. Modern disc mod- 
els, such as those bv iHubenv et all (|2000l . I2001I ). may not 
perfectly explain the shape of quasar big blue bumps (e.g., 
IShang et al.fl2005h . but can provide insight into a number 
of assumptions that go into making bolometric corrections. 
Below we examine a number of features of these state-of-the- 
art disc models and what they tell us about this problem. 

This paper is organized as follows. Section [2] describes 
the model for the continuum spectrum emitted by thin ac- 
cretion discs in type-1 quasars that we use; Section [2]2] ana- 
lyzes the dependence of the disc emission on the viewing an- 
gle and Section [2. 31 outlines how we calculate the bolometric 
luminosity emitted by the disc and the bolometric correc- 
tions. Section [3] contains the theoretical results, with i]3.1l 
describing the dependence of the bolometric correction on 
different wavelengths in the optical-UV band, ^3.2l the corre- 
sponding dependence on the mass accretion rate and black 
hole mass, and i)3.3| giving an account of the systematic error 
obtained in the bolometric correction when quasars are as- 
sumed to emit isotropically. Section |4] discusses the broader 
implications of our work and its possible applications to im- 
prove bolometric corrections using quasar observations. Fi- 
nally, Section [5] summarizes our results. 



2 MODERN THIN ACCRETION DISC 
MODELS 

We intend to perform a self-consistent theoretical calcula- 
tion of the bolometric luminosities of quasars, based on their 
ultraviolet-to-optical continuum spectrum (i.e. the big blue 
bump) and correlate it with the emission properties at differ- 
ent wavelengths. In order to do this, we need detailed model 
spectra of accretion discs, appropriate for quasars. 

The underlying assumption behind the use of the big 
blue bump as a proxy of the bolometric luminosity is that 
the UV-optical band is where most of the energy in the 
quasars is released, whereas emission in other bands is either 
a product of reprocessing of the UV light (e.g. IR) or are 
not en ergetically significant on their own (e.g. radio) (e.g., 
iKoratkar fc Blaes.1999. '). 

Over the years, many models for the continuum spec- 
tra of thin accretion discs have bee n formulated with 
varying degrees of s ophist ication (e.g., iNovikov fc Thornd 
19731: ISun fc Malkan 1989; Laor fc Netzcr 1989: 'Laoij|l990l: 



Sincell fc Krolik 1998 : Hubenv et al.l [2OOO, 200 1]). For our 



urpo ses, we adopted the detailed model of iHubenv et al.l 
2000l ). which consists of a time-steady, geometrically thin, 
optically thick accretion disc which includes general rela- 
tivistic and non-LTE effects. Comptonization is neglected 
since it has little effect for the range of black hole masses 
that we explored l|Hubenv et al.ll2001^ . The Hubeny et al. 
model calculates self-consistently the vertical structure of 
the disc together w ith the radiation field. It was shown by 
IShang et al.l l|2005l ) to reproduce some observed properties 
of the UV spectr a of quasars and is also supported by near- 
IR observations l|Kishimoto et al]|2008l ). These models may 
not be perfect, but they are currently the best theory has 
to offer. 



2.1 Parameter space explored 

iHubenv etall (|2000l ) constructed a grid of models for a wide 
range of values of black hole mass M and mass accretion 
rate M , for two values of the viscosity parameter a (0.01 
and 0.1) and for the two extreme values of the black hole 
spin a corresponding to a maximally spinning Kerr black 
hole (ajM = 0.998) and a nonrotating Schwarzschild hole 
(ajM — 0). The black hole masses explored range from 
1.25 X 10® M0 to 3.2 X 10^°Mq ; the accretion rates range from 
2.4 X lO"** Mq yr~^ to 64 A/q yr~^, or expressed in terms of 
the Eddington accretion rate, between 1.1 x 10~®MEdd and 
0.7A/Edd where the Eddington accretion rate is defined as 
MEdd = 22M/{10^ Mq) M0yr"\ The maximal luminosity 
attained by the models correspond to L/LEdd ~ 0.3. Figure 
[T] shows the M vs. M parameters of the Hubeny et al. grid 
of models for the Kerr case. Hubeny et al. also explored the 
effect of the inclination of the disc with respect to the line of 
sight on the observed spectra, exploring angles ranging from 
nearly face-on {0 ~ 0°) to nearly edge-on discs {9 ~ 90°). 

In our exploration of the parameter space of the accre- 
tion disc model, we restrict ourselves to non-LTE models and 
to accretion rates expected to be appropriate for quasars. 
It is generally accepted that depending on amount of mass 
supplied to the central black hole, different accretion flow 
states are possible. For instance, it is thought that in the 
regime of accretion rates appropriate for bright AGNs and 
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Figure 1. M vs. M parameters of the Hubeny et al. grid of 
models for tlie Kerr case. The solid line illustrates the relation 
m = 0.01 and the dashed line corresponds to m = 1. 



quasars, m > 0.01, the flow is in the thin disc state, while 
for m < 0.01 - the regime appropriate for low-luminosity 
AGNs such as LINERs - the flow would be in a radia- 
tively inefficient or advection -dominated state (ADAF; e.g., 
iNaravan fc McClintocklboOSi ). where we deflne the dimen- 
sionless accretion rate m = M/AfEdd. In our analysis we 
consider only accretion rates with m ^ 0.01. In Figure[Tl we 
illustrate the relation m = 0.01 as the solid line and m = 1 
as the dashed line. 



2.2 The impact of the inclination angle 

Optically thick, geometrically thin accretion discs emit light 
anisotropically, so when calculating their total amount of 
power radiated away we must take that aspect into account. 
For a Newtonian disc, the dependency of the specific flux 
on the inclination of the disc with respect to the line o f 
sight e is simply F„ oc cosg (iFrank. King fc Raind |2002| ). 
The disc model of iHubenv et al.l (|2000D takes into account 
general relativistic effects such as beaming, aberration and 
light bending as well as limb darkening, which introduce 
complicated angle-dependencies in the observed spectrum. 

To illustrate this point, Figure [5] shows the spectra of 
models which differ only with respect to the viewing an- 
gle, with all other parameters held flxed. The specific lumi- 
nosity is calculated as the one an observer along a partic- 
ular viewing angle would see if the source were isotropic, 
i.e. we use the equation — A-Kd^Fv where d is the dis- 
tance. Two sets of models are displayed in that figure, with 
parameters appropriate for quasars - M = K^Mq, M = 
0.25 Mq yr"^ - and bright Seyferts - M = 1.25 x 10® M©, 
M = 0.03 Mq yr"^ - both sets with a = 0.1 and m « 0.01. 
The viewing angles in each set of spectra range from the 
edge-on (lowermost spectrum in each case) to the face-on 
case (uppermost spectrum), with models computed for the 
following values off: 89.4°, 78.5°, 66.4°, 60°, 53°, 36.9° and 
8.1°. Figure [2] illustrates that the viewing angle has a pro- 
nounced impact on the spectra of the thin disc. 

A more insightful way of exhibiting the viewing angle 
dependency of the models can be obtained by plotting the 
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Figure 2. Theoret ical spectra computed with the model of 
iHubenv et"al] ||2000|'| with parameters typical of quasar accretion 
discs {M = IO^Mq, M = 0.25 Mq yr"!; solid line) and bright 
Seyferts (M = 1.25 X 10* Mq, M = 0.03 Mq yr'^; dashed line). 
Each set of models differ only with respect to the viewing angle, 
where the lowermost spectrum in each case corresponds to the 
nearly edge-on angle and the uppermost spectrum represents the 
nearly face-on orientation (see text for the values of 6). 



integrated luminosity emitted by the disc as a function of 
viewing angle. We compute the integrated luminosity as 



Lv dv, 



(1) 



where vq and v\ are respectively the lowest and highest fre- 
quencies at which the thin disc radiates. Li will obviously 
have the same ^-dependency as F^. 

Figure [3] shows how Li depends on 6 for the two fidu- 
cial sets of model spectra mentioned previously (and plot- 
ted in Figure [2}, compared to the Newtonian case in which 
Fv cc cos 9. Clearly, the relativistic model of IHubenv et al.l 
(2000) show significant departures from the simple Newto- 
nian accretion model, especially at large inclination angles. 
Such anisotropic emission must be taken into account when 
estimating the total luminosity of the thin disc and therefore 
the bolometric luminosity of quasars. 



2.3 Theoretical calculation of the bolometric 
correction 

The total luminosity radiated by the disc is independent of 
the viewing angle and is calculated from Li as 



1 

JO 

■k/2 



Li{e) sine de d(f> 



(61) sin 61 d6, 



(2) 



since the grid of models produced by IHubenv et al.l (|200G| ) 
was calculated ranging up to the angle 6 « 90°. 

As mentioned before, a standard practice to estimate 
the bolometric luminosity from quasar SEDs is to use the 
optical-UV portion of the spectrum on the assumption that 
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Figure 3. Dependency of the integrated luminosity Li of the thin 
accretion disc on the viewing angle d. The solid line corresponds to 
the "quasar" fiducial parameters used in Figure [2] and the dashed 
line represents the "Seyfert" fiducial parameters. The dotted line 
illustrates the 6-dependency of the simple Newtonian thin disc 
model for which Fy oc cos 6. Note the log scale. 



it is the energetically dominant component of the radiative 
output and power s the emission at other wavelengths (e.g., 
iNetzer et aLll2007l ). We follow this principle here to estimate 
Lboi and simply assume 

Lhoi = L, (3) 

where L is calculated using Equation [2l 

This formulation for computing Lboi has the virtue that 
it accounts for the anisotropic nature of the emission of the 
thin accretion disc, including all the relevant effects. When 
computing the luminosity of the thin disc, many studies have 
assumed that the disc emits isotropically which introduces 
errors. Figure |4] shows the fractional error incurred when the 
bolometric luminosity is estimated assuming isotropic emis- 
sion, which corresponds to using the equation Lboi = Li as 
opposed to equation O for both the quasar and the Seyfert 
fiducial models mentioned before. The fractional error plot- 
ted corresponds to [Li —L)/L. We can see in this figure that 
errors in the estimate of Lboi as high as « 90% for the face- 
on case can occur when the thin disc is wrongly assumed to 
radiate isotropically. 



3 RESULTS 

In this section, we explore the different results we are able 
to obtain from the grid of models for quasar spectra de- 
scribed in !j2l We are concerned particularly about the re- 
sulting bolometric luminosities and their dependence on the 
properties of the central engine such as black hole mass, 
mass accretion rate and viewing angle. We will also explore 
the relation of Lboi to the monochromatic luminosity at dif- 
ferent wavelengths in the optical- UV band. 
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Figure 4. Fractional error {Li —L)/L incurred when the bolomet- 
ric luminosity is estimated assuming that the accretion disc is an 
isotropic radiator, as a function of viewing angle, for the quasar 
and the Seyfert fiducial models (parameters listed in ^2.21 1. 



3.1 Dependence of bolometric correction on 
wavelength 

One quantity that is of particular interest is the "bolometric 
correction" = Lhoi/{^Lx), which provides a way of esti- 
mating the bolometric luminosity knowing the monochro- 
matic luminosity at a certain wavelength A in the optical- 
UV band. According to the unification scenario of AGNs, 
the optical-UV emis sion would only be directly obse rvable in 
type-1 objects (e.g.. lBartheilll989l , lAntonuccj|l993l '). There- 
fore we restrict our study of the QSO synthetic spectra to 
this class of objects. We do so by restricting the viewing 
angles that we consider to 9 ^ 60°, which we call "type- 
1" models. We note that the opening angle may differ from 
this value, and that it itself may be dependent on luminosity 
(e.g., receding torus models, e.g., Lawrence 1991 ; SimpsonI 
boOSh or other parameters. 

It is worthwhile to analyse the distribution of values of 
("a for different regions of the parameter space, for different 
values of A. We calculate the values of for three wave- 
lengths: 1450 A, 3000 A and 5100 A. Flux measurements at 
these wavelengths are often adopted in observa tional stud- 
ies to make estimates of b olometric luminositie s (|Elvis et al.l 
199|; iKaspi et al] l2000l : iRichards et~al] l2006h . The wave- 
lengths 1450 A, 3000 A and 5100 A correspond to spectral 
regions in the optical for high-redshift, intermediate-redshift 
and low-redshift quasars respectively, and at these wave- 
lengths the emission lines are not present or are relatively 
weak. 

Figure \S\ shows the histograms of the distribution of 
values of for non-LTE models corresponding to "type-1" 
models with m ^ 0.01, i.e. "type-1 quasars". We can see in 
Figure [S] that as the wavelength increases, the distribution 
of values of broadens. The peak of the distribution moves 
towards higher values as A increases, with the medians of the 
bolometric correction for each wavelength being (Ci45o) ~ 3, 
(Csooo) = 5.9 and (Csioo) = 7.6. 

Figure |6] shows the relation between Lboi and ALa with 
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Figure 6. Relation between Lbol a^nd -^^A at A = 5100 A, for 
the "type-1" quasar model spectra. The dotted and dashed lines 
correspond to models of the form Lbol = CsiOO^^Ai while the 
solid line show the best-fitted power-law to the data using the 
method. 

A — 5100 A, for the type-1 quasar models with m 0.01. 
The dotted and dashed lines in Figure [6] correspond to 
models in which the bolometric luminosity is simply re- 
lated to the monochromatic luminosity at 5100 A through 
^boi = CsiooAZ/A, with the value of ("sioo being respectively 9 
and 13 for the dotted an d dashed lines. The se va lues of Csioo 
are the ones derived bv lKaspi etatl (|2000l ) andlM vis et al.l 
respectively. The mean deviation (T(logLboi) of the 
data about these lines is 0.47 (Csioo = 9x model) and 0.49 
(Csioo = 13x model). 

We parametrize a power-law fit in the log-log space por- 
trayed in Figure [6] as logLboi ~ A + Blog{XLx) and fit it 
to the data using the statistic, obtaining the parameters 
A = 11.7±0.93 and B = 0.76±0.02 (solid line), with a mean 
deviation about the best-fit model of 0.38. This result sug- 
gests that the non-linear model of the form Lboi oc (ALa)"'^^ 
is favoured to explain the data for the sample of simulated 
type-1 quasar spectra for A — 5100 A. 



Figure 7. Same as Figure Elfor A = 3000 A. 

Figure [7] shows the Lboi vs. XLx for type-1 quasar mod- 
els at A = 3000 A, following the same notation as Figure 
El A fit to the 3000 A data using the method results in 
A = 9.24 ±0.77 and B = 0.81 ±0.02 (dot-dashed line), with 
a mean deviation about the best-fit model of 0.31. For com- 
parison, the corresponding mean deviations about the linear 
models with (^3000 ~ 9 and ^3000 = 13 are respectively 0.38 
and 0.46. 

Figure [8] shows the Lboi vs. XLx for type-1 quasar mod- 
els at A = 1450 A, again following the notation of Figure 
[H A fit to the 1450 A data using the x^ method results in 
A = 6.7 ± 0.69 and B = 0.87 ± 0.02, with a mean devia- 
tion about the best-fit model of 0.26. For comparison, the 
mean deviations about the linear models with 1^1450 ~ 9 and 
C1450 = 13 are respectively 0.51 and 0.64. 



3.2 Dependence of bolometric correction on the 
accretion rate and black hole mass 

In this section we analyze the possible dependence of the 
bolometric correction on the Eddington ratio Lboi/LEdd, 
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Figure 10. Same as Figure O but here the filled circles cor- 
respond to M < 10® Mq and the open circles correspond to 
M ^ 10® Mq. 
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Figure 9. Relation between Lbol and \L\ at A = 5100 A, for 
the "type-1" quasar model spectra. The filled circles correspond to 
^bol/^Edd < 0-1 S'lid the open circles correspond to ibol/iEdd 
0.1. The solid and dot-dashed lines correspond to a broken power- 
law fit and a simple power-law, respectively. 



which traces the accretion rate, and the black hole mass 
for the "type-1" quasar model spectra. 

Figure |§] shows the relation between Lboi and \L\ 
at A = 5100 A, where the filled circles correspond to 
iboi/ifidd < 0.1 and the open circles correspond to 
iboi/^Edd ^ 0.1. There is a hint in this plot that the slope 
of the dependence of Lboi on \L\ changes slightly as the 
Eddington ratio increases, with a break at log(ALA) ~ 45.5 
corresponding to the chosen "transition value" I/boi/i'Edd = 
0.1. The solid line in Figure [9] displays the broken power- 
law best-fitted to the data using the Marquardt-Levenberg 
least-squares method, with the break at a; = 45.5. The mean 
deviation of the data about this best-fit model is 0.35. 

Figure [10] is similar to the previous one, but now the 
filled circles correspond to models for which M < lif Mq 
while the open circles represent models with M ^ 10® Mq. 
We find a slight hint in this figure of a change in the slope of 
the power-law fit as the black hole mass increases, similarly 



to what was observed in the previous plot. The solid line 
shows the same broken power-law fit shown in Figure [§1 

There is no strong reason to abandon a simple power- 
law model in favour of more complicated models such 
as a broken power-law, since the latter reproduces only 
marginally better the results compared to a simple power- 
law. 



3.3 Beyond the isotropy assumption when 
computing quasar luminosities 

When converting quasar nuclear line-of-sight fiuxes to 
luminosities or when computing bolometric corrections 
and luminosities, the usual assumption is that quasars 
emit isotropically (e .g., lElvis et al.lll994l : IShang et aLll2005l : 
[Richards et al.ir2006l ). We know that this is not the 
we showed in 8^2.31 (c.f. Figure |4)|. Using the spectral models 
for quasar accretion discs, we can have a theoretical per- 
spective on how good (or bad) is the isotropy assumption 
on average. 

In order to address this issue, we calculated the average 
error implied by the isotropy assumption when estimating 
the bolometric luminosities of accretion discs from line-of- 
sight fiuxes measured in the optical-UV continuum spec- 
trum. We did this by averaging all the type-1 quasar models 
at specific viewing angles. Figure [Tl] shows the correspond- 
ing mean fractional error calculated in this way, defined as 



Mean fractional error - 



Lj — Lb 

Lhn] 



(4) 



(see equations [T] and 131) as a function of the viewing angle 0. 
We show in this figure the fractional error that result when 
two different accretion disc models are used: the general rel- 
ativistic model (toget her with the la un certainty for the 
relativistic model) of Hubenv et al.l ||2000|) and for compar - 
ison a simple Newtonian disc I Frank. King fc Rainell2002l ). 
For the Newtonian model, the bolometric luminosity is re- 
lated to the integrated luminosity assuming isotropy Li as 
Lboi ~ 1/(2 cos8)Li, and the fractional error is given simply 
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Figure 11. Average fractional error incurred when the bolomet- 
ric luminosity is estimated assuming that the accretion disc is an 
isotropic radiator, as a function of viewing angle. 



by 2 cos — 1. For the general relativistic model, we calculate 
the mean error numerically. 

For each accretion disc model there is a specific line-of- 
sight 9i at which Lboi = Li, i.e. at these viewing angles the 
mean error is zero and there is no error incurred in using the 
isotropy assumption. For the Newtonian case, this angle is 
simply 9i = 60°. For the relativistic case. Figure [TT] tells us 
that 9i ~ 66°, with the uncertainty in the mean error being 
lowest at S « 55°. 



4 DISCUSSION 

We note that the three main underlying assumptions be- 
hind our work are (1) that the models for the spectra of 
thin accretion disc correctly explain the observed optical- 
UV continuum spectra of quasars, (2) the parameter space 
for the accretion disc models that we explored closely resem- 
ble that of real quasars, and (3) the X-ray contribution to 
the bolometric luminosity is small compared to that of the 
big blue bump. The resulting scatter of (T(log Lboi) ~ 0.9 dex 
in the model data (c.f. Figure O arises from the variation 
in the different parameters of the central engine: accretion 
rate, black hole mass and spin, viewing angles and viscos- 
ity. The fact that the theoretical scatter is similar to the 
scatter measured in observed quasar data (|Elvis et al.lll994 
IShang et al .120051 : [Richards' et al ]|2006l ) suggests that our as- 
sumptions 1 and 2 are reasonable and offers an explanation 
to the observed scatter. 

There is accumulating evidence that less than ~ 3% 
of the radiative output of quasars is r e leased in X-rays 
jElvis et al] Il994l: rv^_s udevan fc FabianI l2007l : IHoI l2008l : 
IVasudevan fc Fabian 2009). This indicates that our third 
assumption is quite plausible. 

Previous results based on observed quasar SEDs find 
simple linear bolometric corrections of the form I/boi oc 
XL\ for different reference wavelength s (|Elvis et al ] 11994 
iKaspi et al]|2000l : iRichards et all 120061 '). Contrary to that, 
we find that non-linear bolometric corrections of the form 



log Lboi = A + Blog{XLx) better reproduces the synthetic 
data resulting from the accretion disc models, with B ^ 0.9 
for the wavelengths in the optical-UV band that we con- 
sidered. Our findings are consistent with the observations 
which s how that the quasa r SED shape changes with lumi- 
nosity l|Steffen et al.ll2006l ). On the other hand, we found 
only weak evidence for evolution of the bolometric correc- 
tion with the Eddington ratio and the black hole mass. 

We carried our calculations of the bolometric correc- 
tions for three fiducial wavelengths, 1450 A, 3000 A and 
5100 A. We find that as the wavelength decreases, the dis- 
tribution of values of L]:,oi/{XLx) becomes narrower (c.f. 
Figure [5]). Furthermore, the relation between Lboi and 
the monochromatic luminosities becomes tighter, in the 
sense that both the uncertainties in the parameters of the 
power-law fits and the dispersion around the best-fit model 
decrease f ^3.1|l . This result simply refiects the fact that 
monochromatic luminosities measured closer to the peak 
of the energy distribution tend to be better tracers of the 
total energy released by the accretion disc. Although it is 
preferable to use luminosities emitted at increasingly smaller 
wavelengths for bolometric corrections, there is one impor- 
tant practical caveat as we advance upwards in energy: the 
effects of dust extinction become progressively more strong, 
which obviously must be considered. 

We can make some recommendations based on the 
models for the spectra of accretion discs that we stud- 
ied, which can be used to estimate bolometric luminosities 
from monochromatic luminosities at specific wavelengths. 
For A — 1450 A, the following bolometric correction formula 
is suggested from the models (see ^STTJ: 



log Lboi = (6.7 ± 0.69) + (0.87 ± 0.02) log L1450 • (5) 
For A = 3000 A, we derived the correction 

log Lboi = (9.24 ± 0.77) + (0.81 ± 0.02) logLgooo- (6) 
Finally, for A = 5100 A, the following correction is indicated: 

log Lboi = (11.7 ±0.93) + (0.76 ±0.02) log L5100. (7) 

Since the bolometric corrections are tighter for A — 1450 A, 
it is preferable to use the measurements at this wavelengths. 
When using these bolometric corrections, the readers should 
be aware of the underlying assumptions of our work de- 
scribed previously. 

We demonstrated that the widespread assumption that 
quasars emit isotropically leads to considerable errors in the 
determination of the actual bolometric luminosities, with 
average fractional errors exceeding 30%. How can these the- 
oretical results be conceivably applied in a real-world case, 
to improve observational estimates of the quasar bolomet- 
ric luminosity? Let us assume that the bolometric luminosity 
Lb°i is available from integrating an observed quasar optical- 
UV spectrum adopting isotropy, but the line-of-sight to the 
quasar is unknown (as is usually the case). If all lines of 
sight to the central engine are in principle available (i.e. no 
dusty torus), then the average viewing angle for a quasar 
is (0) « 57° . Figure [TT] tells us that the mean systematic 
fractional error is 0.1 and the actual bolometric luminosity 
would be systematically smaller than the isotropic estimate 
by a factor of « 0.9. 

Of course, we must consider the presence of a dusty 
torus blocking our view to a large part of the accretion disc. 
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If we consider a dusty structure such that only viewing an- 
gles 6 < 60° are accessible, then (6) ~ 39°. In this case, 
the systematic average fractional error resulting from the 
isotropic estimate is 0.24, such that Lboi ~ 0.81/^°!. A some- 
what larger dusty torus that restricts the lines of sight to 
6 < 45° would correspond to {9} « 30°; the corresponding 
average fractional error would be 0.28 with Lboi ~ 0.781/boi- 
It is important to keep in mind these systematic errors 
when estimating quasar bolometric luminosities assuming 

isotropy. 

In the cas e of receding torus models (|Lawrencelll99ll : 
ISimpsonI [20051 ). which have larger opening angles with in- 
creasing quasar luminosity, edge-on discs become more com- 
monly seen in the most luminous systems. In this case, bolo- 
metric corrections would be expected to be a function of 
luminosity. Additionally, the isotropy assumption would be 
less in error for a larger fraction of luminous quasars, on 
average. 

Our results indicate that determining torus opening an- 
gles and their distribution as a function of luminosity would 
help improve bolometric corrections. Also finding ways of 



determinin g disc inclination angles, p erhaps through radio 
properties (Wills fc Brothertonl l l995l ) or other approaches 
l|Down et al.ll2010f ) is also indicated as a way of significantly 
improving bolometric corrections in individual objects. 



5 SUMMARY 

We have studied state-of-the-art models for the continuum 
spectrum of thin accretion discs, exploring the theoretical 
optical-to-UV bolometric corrections for a range of mass ac- 
cretion rates, black hole masses and viewing angles appro- 
priate for the case of type-1 quasars. We find a number of 
results of direct relevance in observational studies aimed at 
the determination of quasar bolometric luminosities. 

Regarding the relation between the bolometric and the 
monochromatic luminosities measured in a given wavelength 
in the optical-to-UV band, we obtain: 

(i) A nonlinear relationship of the form Lboi oc (ALa)°°"''* 
is favoured by the models, with the deviation from linearity 
increasing as the wavelength decreases. 

(ii) The bolometric corrections become tighter with de- 
creasing wavelength, i.e. as we approach the peak of the 
quasar SED (although reddening affects shorter wavelengths 
more strongly). 

(iii) We list fitting formulas for the bolometric correction 
corresponding to the wavelengths A = 1450 A, 3000 A, and 
5100 A, which can be used to improve observational esti- 
mates. 

Quasars are not isotropic emitters, and making this as- 
sumption may lead to severe systematic errors in the deter- 
mination of Z/boi using the method of integrating the "big 
blue bump" spectrum. For a typical line-of-sight of ~ 30° or 
less to the quasar central engine, we estimate that there is 
a systematic error of ~ 30% or higher implied if Lboi is esti- 
mated using the isotropy assumption. We find that a more 
careful calculation should be Lboi ~ O.SLboi) where Lboi is 
the bolometric luminosity estimated from integrating the 
optical-to-UV spectrum assuming isotropy. 

With the results of this work, a better assessment of 



the radiative output of quasars based on observations will 
be made possible. 



ACKNOWLEDGMENTS 

We thank Omer Blaes for making the library of thin disk 
spectral models available to the community. RSN is grateful 
to the hospitality of Department of Physics & Astronomy at 
the University of Wyoming, were part of this work was car- 
ried out, and acknowledges the financial support of CNPq. 
MB is grateful to the hospitality of Instituto de Fi'sica, Uni- 
versidade Federal do Rio Grande do Sul in Porto Alegre, 
Brazil, and acknowledges support from NASA through grant 
NNG05GE84G. 



REFERENCES 

Antonucci, R. 1993, ARA&A, 31, 473 
Barthel, P. D., 1989, ApJ, 336, 606 

Begelman, M. C. 2004, in Coevolution of black Holes and 

Galaxies, Carnegie Observatories Astrophysics Series, ed. 

L. C. Ho, p. 374 
Crenshaw, D. M., Kraemer, S. B., George, I. M., 2003, 

ARA&A, 41, 117 
Di Matteo, T., Springel, V., & Hernquist, L. 2005, Nature, 

433, 604 

Down, E. J., Rawlings, S., Sivia, D. S., Baker, J. G., 2010, 

MNRAS, 401, 633 
Elvis, M., et al. 1994, ApJS, 95, 1 

Ferrarese, L., & Ford, H. 2005, Space Sci. Rev., 116, 523 
Frank, J., King, A. R., & Raine, D. J. 2002, in Accretion 

Power in Astrophysics (3d ed.; Cambridge: Cambridge 

Univ. Press) 
Ho, L. G. 2008, AAR&A, 46, 475 

Hubeny, I., Agol, E., Blaes, O., & Krolik, J. H. 2000, ApJ, 
533, 710 

Hubeny, I., Blaes, O., Krolik, J. H., & Agol, E. 2001, ApJ, 
559, 680 

Kaspi, S., Smith, P. S., Netzer, H., Maoz, D., Jannuzi, B. 

T., & Giveon, U. 2000, ApJ, 533, 631 
Kishimoto, M., Antonucci, R., & Blaes, O. 2005, MNRAS, 

364, 640 

Kishimoto, M., et al. 2008, Nature, 454, 492 

Koratkar, A., & Blaes, O. 1999, PASP, 111, 1 

Laor, A., & Netzer, H. 1989, MNRAS, 238, 897 

Laor, A. 1990, MNRAS, 246, 369 

Lawrence, A., 1991, MNRAS, 252, 586 

Narayan, R., & McChntock, J. E. 2008, New Astronomy 

Review, 51, 733 
Netzer, H., et al. 2007, ApJ, 666, 806 

Novikov, I. D., & Thorne, K. S. 1973, in Black Holes, ed. C. 
De Witt & B. S. De Witt (New York: Gordon & Breach), 
344 

Richards, G., T., et al. 2006, ApJ, 166, 470 

Shang, Z., et al. 2005, ApJ, 619, 41 

Simpson, C, 2005, MNRAS, 360, 565 

SinceU, M. W., & Krolik, J. H. 1998, ApJ, 496, 737 

Steffen, A. T., et al. 2006, AJ, 131, 2826 

Sun, W.-H., & Malkan, M. A. 1989, ApJ, 346, 68 



Quasar bolometric corrections 9 

Vasudevan, R. V., & Fabian, A. C. 2007, MNRAS, 381, 
1235 

Vasudevan, R. V., & Fabian, A. C. 2009, MNRAS, 392, 
1124 

Vestergaard, M., 2009, Spring Symposium on Black 

Holes. Cambridge Univ. Press, Cambridge, in press 

l|ajXi^ :0904.2615 ) 
Wills, B. J., Brotherton, M. S., 1995, ApJ, 448, L81 



This paper has been typeset from a T^rjX/ I^T^jX file prepared 
by the author. 



